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Abstract 

We test the ability of large scale velocity helds inferred from geomagnetic secular 
variation data to produce the global magnetic field of the Earth. Our kinematic dynamo 
calculations use quasi-geostrophic (QG) flows inverted from geomagnetic field models 
which, as such, incorporate flow structures that are Earth-like and may be important 
for the geodynamo. Furthermore, the QG hypothesis allows straightforward prolongation 
of the flow from the core surface to the bulk. As expected from previous studies, we 
check that a simple QG flow is not able to sustain the magnetic field against ohmic 
decay. Additional complexity is then introduced in the flow, inspired by the action of the 
Lorentz force. Indeed, on centenial time-scales, the Lorentz force can balance the Coriolis 
force and strict quasi-geostrophy may not be the best ansatz. When our columnar flow 
is modihed to account for the action of the Lorentz force, magnetic held is generated 
for Elsasser numbers larger than 0.25 and magnetic Reynolds numbers larger than 100. 
This suggests that our large scale how captures the relevant features for the generation 
of the Earth’s magnetic held and that the invisible small scale how may not be directly 
involved in this process. Near the threshold, the resulting magnetic held is dominated 
by an axial dipole, with some reversed hux patches. Time-dependence is also considered, 
derived from principal component analysis applied to the inverted hows. We hnd that 
time periods from 120 to 50 years do not affect the mean growth rate of the kinematic 
dynamos. Finally we notice the footprint of the inner-core in the magnetic held generated 
deep in the bulk of the shell, although we did not include one in our computations. 


1 Introduction 

The main Earth’s magnetic field and its temporal variations are generated by the motions of 
liquid metal in the core. Provided some assumptions are made, it is possible to infer the large 
scales of the flow at the top of the Earth’s core, from observations of the geomagnetic field and 
its variations with time, which are commonly referred to as Secular Variation (SV). Because 
the crustal magnetic field dominates at small spatial scales, the core field is known only for 
the largest scales (up to spherical harmonic degree 13). Similarly, because of time varying 
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currents in the magnetosphere and in the ionosphere, the SV produced by core processes can 
be isolated only up to harmonic degree 12 to 14 [e.g. Olsen et al., 2014 . Unfortunately, this 


inherently also limits to relatively large scales the flow we can reconstruct at the top of the 
core. 


Direct numerical simulations of the geodynamo, pioneered by Glatzmaier and Roberts 


1995 


can be tuned to produce magnetic fields that resemble closely the one of the Earth [e.g. 

However, the mechanism 


Christensen et ah, 2010| Christensen, 2011, Aubert et ah, 2013 


by which the magnetic field is actually generated in the Earth’s core is not clear. Indeed, 
the parameter range in which those simulations operate is arguably very far from the one 
expected in the Earth, and when extrapolating them to realistic Earth parameters, we obtain 
a picture where the convective flow is important down to small scales of 0.1 to 10 km wide. 
Is this hidden small scale flow (both from numerics and from inversions) important for the 
generation of the magnetic field? Or is the large scale flow alone responsible for the induction? 
To help answer these difficult questions, the modest goal of this paper is to test the capability 
for dynamo action of the large scale flows inferred from geomagnetic field models. 

The importance of the Coriolis force in the core arguably leads to flows that are mainly 
invariant along the rotation axis, which are referred to as quasi-geostrophic (QC) or columnar 
flows. These columnar flows are thought to be relevant at large spatial scales and short time- 
Cillet et al., 2011|. At smaller spatial scales or more importantly at longer 


scales Jault,2008 


time scales, the columns are expected to wither: at such scales the Lorentz and buoyancy 
forces should arguably be both taken into account. The effect of the Lorentz force is of 
particular interest as Sreenivasan and Jones 2011| have shown that the induced magnetic 
pumping significantly enhances the dynamo action of the flow. 

In this study, we use a columnar flow assumption to infer the flow at the core surface. 
The reconstruction of the flow inside the whole core is thus straightforward. We then use this 
reconstructed flow in a kinematic dynamo code to explore its dynamo capability. Despite being 
a much simpler approach than solving the full dynamical system, solving only the induction 


equation has already proven to be useful [e.g. Cubbins, 2008 Tobias and Cattaneo, 2013 


Cabanes et al., 2014 


Although care must be taken when applying anti-dynamo theorems 
see e.g. Bachtiar et ah, 2006, Li et ah, 2010] , the quasi-two dimensionality 

Despite this 


to bounded flows 

of our columnar flows does not a priori help dynamo action jjones 


2008 


fact, columnar flows have already proven to be capable of dynamo action in the presence 
of large-scale zonal shears [Schaeffer and Cardin 2006 Cuervilly, 2010 . In addition to the 
large-scale shears, both studies stress the importance of time-variability of the flow for the 
dynamo action. As for inverted core flows, they never exhibit dominant large scale shears. 
To try to boost the dynamo efficiency of columnar flows, we introduce a time-dependence 
matching inferred core flows and perturbate the flow with the magnetic pumping described 
by [Sreenivasan and Jones 2011 . As we will show, only the latter actually leads to dynamo 


action. 

Even though our flow explains an important part of the observed SV, we do not expect 
the growing dynamo field to match any particular feature of the secular variation as, e.g., the 
dipole moment decay rate. The reasons are inherent to the kinematic dynamo problem and 
are detailed in section [5l 

The point we wish to highlight in the present study is that our flow model, obtained 
through a geomagnetic field model inversion, and including a physically relevant magnetic 
pumping that does not change the surface flow, is efficient to maintain the magnetic field of 
the Earth, without the need for contributions from smaller scales. 
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The paper is organized as follow: in the next section the method for obtaining the core 
flows is described and discussed. Then we describe the procedure we follow to compute 
kinematic dynamos from the core flows, and introduce the relevant parameters. The results 
are first described in section and then discussed in section We end the paper with 
concluding remarks. 


2 Inverted core flows 

2.1 Columnar flows 

Columnar flow can be seen as the outcome of competing mechanisms that propagate infor¬ 
mation inside the core: Alfven waves due to the magnetic tension in the medium and inertial 
waves due to the restoring effect of Coriolis force. On short time-scales and large length-scales, 
these latter waves are quicker in propagating perturbations axially (along planetary vorticity 
lines) inside the liquid core, as expressed by a small Lehnert number, i.e. the ratio between 


2015 


Alfven wave and inertial wave speeds Jault, 2008, Gillet et ah, 2011, Nataf and Schaeffer 


Schaeffer and Pais 2011 pointed out that equatorially anti-symmetric (AS) features can 


be present inside a spherical rotating container as a result of some anti-symmetric excitation. 
When looking for this symmetry in core flows inverted from geomagnetic field models, they 


further noticed that the AS component tends to prevail for small scale flows [see also Gillet 


et ah, 2011] , which can be understood in terms of a lengthscale-dependent Lehnert number, 
Xi = B/iQy/pJIo, larger for small-scale structures. 

In a recent work. Pais et al. [2015 applied Principal Component Analysis tools to the 
streamfunction ^ of the QG flow inverted from geomagnetic field models. The ‘data’, which 
consisted in the values of streamfunction ^ on a regular spatial grid at the core surface, were 
decomposed into a mean flow plus a linear combination of a small number of spatial patterns 
multiplied by time-varying coefficients. This approach is particularly useful in the present 
context, by providing the means to describe the time-varying flow derived from geomagnetic 
field model gufml over a 150 year time period with a small number of parameters (see section 


3.3.1). 


Since allowing for an AS component in inverted core flows improved variation of Length- 

argued in favor of considering both 


of-Day (ALOD) estimations, Schaeffer and Pais 2011 
flow symmetries in the inversion. 


2.2 Obtaining the core flow 

The geomagnetic field model gufml has been inverted for a large-scale columnar flow and 
an AS component, providing for a flow solution for each epoch in the 1840-1990 period (one 
independent flow snapshot each year). The inversion acconnts for a separation of QG and AS 
flows as in Schaeffer and Pals] |2011j . 

The columnar flow verifies the following kinematic condition at the core-mantle boundary 
(GMB): 

= ( 1 ) 

with Vc the Earth’s core radius, whereas the AS component satisfies 


Vh • u =0 


( 2 ) 
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Figure 1: On the top row, the mean flow using a Hammer projection centered at the 180° 
meridian (left) and an ortographic projection as seen from the North pole (right). On the 
bottom row, the three spatial structures or Empirical Orthogonal Functions characterizing the 
first three modes from the PC analysis. The EOF’s are normalised and the scale is arbitrary. 
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These are the surface constraints of a flow that in the bulk of the core has a QG component, 
equatorially symmetric (upper index ‘+’)) an AS component (upper index 

In Schaeffer and Pais |2011| , the AS component u had no particular kinematical con¬ 
straint imposed. By using ([^ instead, the AS surface flow component can be completely 
retrieved from a toroidal scalar T, in the same way as the QG surface flow component can 
be retrieved from a streamfunction scalar ^ [see e.g. 


Pais and Jault 2008, Pais et ah, 2015 


u+ = 


1 


AC(0,0)r + 


sm ( 


cos 6 

u~ = 'Vh/\T{ 0,4>)r 


Tr COS^ 9 




(3) 

(4) 


where ^ and T are given in rad/yr. In this study, where we test the ability of QG flows to 
increase the energy of magnetic field modes that have a similar morphology to the Earth’s 
field, we extracted u"*" from the inverted flows and did not further consider the AS component 

u“. 

Besides the two conditions Q and ([^ imposed on the surface flow through quadratic forms 
on the flow coefficients multiplied by relatively high regularization parameters, two further reg¬ 
ularizations are used: a penalization of the azimuthal gradients, /( 7 j\^g[(l/sin0) {du/d(l))^]dS 
as in Schaeffer and Pais [2011 , supporting the observation of structures developed preferably 
along parallels in natural rotating flow systems; a penalization of radial vorticity and hori¬ 
zontal divergence, H A u)^ -|- (Vh ■ u)^]d5 , corresponding to the norm used in 

Gillet et al. |2009| and favoring a large scale flow. 


Following Pais et al. [2015 , the streamfunction ^ obtained for the QG flow is analyzed 


into a linear combination of a small number of spatial patterns multiplied by time-varying 
coefficients, using Principal Component Analysis (PGA). Five modes are retained, which 
explain more than 95% of the time variability of the flow. The reduced streamfunction model 
reconstructed from these main modes is given by 


C(rc, 9, (j), t) = ^°(rc, 9,(j)) + '^ PCk{t) CEOpirc, 9, (p) 


(5) 


k=l 


where is the mean flow, obtained by averaging the flow coefficients over the time-period 
1840 - 1990, Empirical Orthogonal Function (EOF) of order k depending on 

spatial coordinates, and PCk{t) is a function of time, the Principal Component (PC) of order 
k. The product PCk{t) ^^Qp{rc,9,(f>) is the A:th PGA mode, and explains a percentage of 
the time variability found in data. The different modes are uncorrelated in time during the 
1840-1990 interval and in space over the whole CMB. Figure shows the flow at the CMB 
captured by and the first three EOFs. Figure shows plots of the five first PCs that enter 
expression [sj Note the resemblance between this flow and those in Schaeffer and Pais 2011| 


and Pais et al. [2015 


In all cases, the mean flow shows strong azimuthal currents at high 
latitudes centered at -145° longitude and at low latitudes centered at 0°. The three large 
vortices at medium/high latitudes described in Pais et al. [2015 are also present. 

The polar anticyclone lying inside the tangent cylinder cannot be conveniently retrieved 
with conventional QG flow regularization 


[e-g- 


Pais and Jault, 2008 


It is nonetheless an 

important feature of core flows which has already been discussed in both observations [e.g. 


Olson and Aurnou, 1999 , and geodynamo simulations [e.g. Aubertf 2005 . 
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3 A kinematic dynamo problem 

3.1 Induction equation 

The evolution of the magnetic held within an electrically conducting huid is given by the 
induction equation 

= V X (v X B) + r/V^B. (6) 

at 

where rj = is the magnetic diffusivity, is the magnetic permeability of empty space, 

and a is the electrical conductivity of the huid that we assume to be homogeneous. Here, v 
denotes the entire three-dimensional how in the bulk, while u is the how at the core surface. 
The key parameter for dynamo action is the magnetic Reynolds number 

Rm = VoTc/'r^ = Vorcfioa 


which compares the magnitude of the induction term to the ohmic dissipation. We use as 
characteristic speed Vq the maximum value of the velocity of the how held v, and as length 
scale the radius of the core. The core hows in this study have Vq ~ 15 km/yr, leading to 
Rm ~ 800 to 3000 for the Earth’s core, depending on the estimated electrical conductivity 


P 


ozzo et al. 


2012 


Dynamo action, which is the spontaneous growth of a magnetic held from the motion of a 
conducting huid, happens when the induction overcomes ohmic dissipation, which occurs for 
Rm > Rmc- Numerical computations are generally needed to determine the critical magnetic 
Reynolds number Rmc, because it depends on the precise details of the how. Efficient hows 

2008| . Note however. 


in spheres have Rmc ~ 10 to 100 [e.g. Dudley and James 1989, Jones 


that not all how helds can ultimately trigger dynamo action, in which case Rmc is not dehned. 
Einding Rmc and the fastest growing magnetic held B when the how v is prescribed is referred 
to as the kinematic dynamo problem. 

2008| and Sarson [2003 


In the context of geomagnetism, Gubbins 


studied the kinematic 

dynamo problem for minimalistic hows that captured some expected features of core hows. 
Here, the prescribed velocity held v and its time evolution are expected to be more Earth-like, 
since they are derived from the core how inverted from geomagnetic held observations as 
described in section [H This will be detailed below. 


3.2 Numerical procedure 


In order to determine Rmc for a given how v, we compute the time-evolution of B as given 
by the induction equation Q, for several values of the magnetic diffusivity rj. 


et al. 

2011 

Monteux et al. 

2012 

Cabanes et al. 

2014 


and has been adapted for this study [see Schaeffer et al 
transform library SHTns 


Schaeffer, 2013 


2015 


which is available as free software. 
It uses the spherical harmonic 
and hnite differences in radius. The integration 
scheme is second order in time with the diffusion terms treated by a Crank-Nicolson scheme, 
whereas an Adams-Bashforth one is used for the induction term. 

Both the prescribed flow v and the magnetic held B occupy the full sphere: we have no 
solid inner-core in our computations to avoid the issues arising when trying to reconstruct a 
columnar how compatible with a solid inner-core. Note that the XSHELLS code has been 
benchmarked to full-sphere solutions [Marti et al., 2014[ and used in a previous study involving 


6 










































































full-spheres Monteux et al., 2012| . The mantle is assumed electrically insulating, so that the 
magnetic field in the liquid sphere matches a potential field outside the conducting domain. 

The magnetic energy is monitored, its growth rate indicating if the flow leads to dynamo 
action. However, when starting from a random magnetic seed, transient decay and growth 
may be observed before the average growth-rate converges to a constant value. The time 
needed to obtain a well converged mean growth-rate is typically a fraction of the magnetic 
diffusion time = r'^/rj. All our simulations have been running at least for one fifth of 
but often one half of and, in a few cases, up to 2T^. 

A typical run is set up with = 300 radial levels, and spherical harmonic expansions 
truncated at degree ^rnax = 149 and order rumax = 85. For the most demanding cases (largest 
Rm) and to check numerical convergence, we pushed resolution up to {Nr, imax,'mmax) = 
(320,250,150). The computations are always fully dealiased using the appropriate numbers 
of grid points in latitudinal {Ng > ^imaxl‘^) and longitudinal (A^,^ > Srrimax) directions. 


3.3 Prescribed flow models 
3.3.1 Time dependence through PC 


It has been shown that the time-dependence can be important for dynamo action |Willis and 


Gubbins, 2004, Schaeffer and Cardin, 2006, Tilgner, 2008| . We thus also use flows varying 
in time in our kinematic dynamo study, as captured by geomagnetic observations. The 
prescribed velocity field u is computed from equation Q , with ^ (t) given by equation ([^ . In 
order to asses if a flow is capable of sustaining a magnetic field against ohmic dissipation, we 
need to integrate the induction equation ([^ for times comparable to the magnetic diffusion 
time. Because the magnetic diffusion time in the Earth is much longer than the time-period 
for which the core flow u is computed, we cannot reconstruct the flow for a time long enough 
for a simulation to run. To overcome this, we fitted each of the first five principal components 
PCk{t) (see eq. with one sine function PCk{t) = Aksui{ujkt + Uk), as represented in figure 
This allows us to compute a flow u at any time, by extrapolating from decade time scales 
at which observations take place to much larger diffusive ones at which dynamo action may 
occur. 

The prescribed bulk flow v(t) is then reconstructed from the symmetric surface flow u'*'(t) 
using the approach we now describe. 


3.3.2 Purely columnar flow 

The simplest way to reconstruct the flow at any depth in the core, is to assume a columnar 
flow. The cylindrical components {s,(j),z) of the bulk flow v are related to the spherical 
components {9,4>) of the surface flow u"*" by 

Vsis,(j)) = ^^^uj{e,^) 

T c 

V^{s,(p) = u^{0,4>) 

SZ 

v^{s,(p,z) = u-^{9,4') (7) 

with s = rcSinO the cylindrical radial coordinate and H{s) = — §2 iialf height of a 

fluid column. Because u"*" derives from a pseudo-streamfunction (see eq. , mass conservation 
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Figure 2: The main five PCs (used to reconstruct the flow using eq. and their fitted sine 
functions, each defined by an amplitude, a time period and a phase shift. 


is ensured and we always have V.v = 0 [see also Amit and Pais, 2013 . Remember also that 
in this study, we ignore the AS component u~. 


3.3.3 Magnetic pumping 

Helicity is another important ingredient for dynamo action as it leads to the so-called alpha- 
effect [e.g. Jones, 2008 whereby poloidal magnetic field is produced from toroidal field. Ekman 


pumping in columnar vortices produces helicity Busse, 1975 , but it vanishes at small enough 


Ekman numbers [e.g. Schaeffer and Cardin, 2006 


On centenial time-scales the influence of the Lorentz force on the flow should be taken 
into account. Indeed, when a magnetic field permeates columnar vortices, [Sreenivasan and 
Jones [2011 have put forward an effect coined magnetic pumping, which produces an axial 
flow in phase with axial vorticity, that significantly enhances the mean helicity of the flow and 
consequent dynamo action. Because helds of dipolar symmetry have a more efficient mag¬ 
netic pumping, they argued that this effect could explain the preference of dipole-dominated 
magnetic fields in their simulations. 

The magnetic pumping is proportional to the local vorticity and to the square of the 
magnetic field. Sreenivasan and Jones [2011] have computed magnetic pumping solutions for 
a simple toroidal held 

s z(77(s) 2 - z2) 


— Bn 


H{s) 


( 8 ) 


of dipolar symmetry. 
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Because in our kinematic dynamo approach the flow is prescribed by the streamfunction 
? (eq. I and[^, we cannot take into account the retroaction of the dynamic magnetic field. 
Instead, we assume the Earth permeated by a simple toroidal field of dipolar symmetry, and 
follow [Sreenivasan and Jones 2011] , As no explicit expression of the magnetic pumping 
is available, we choose to model the corresponding axial velocity with the following 


parametrization that closely mimics the magnetic pumping computed by Sreenivasan and 
Jones poTI I (see their equation 3.9 and their figure lb). Using a Fourier decomposition 

we prescribe, for all azimuthal wavenumber m: 


v^p = A Vq f{z/H{s)) b{s) rr?irn{s) 


(9) 


where f{x) = —^x{l — x)^(l + x)^ captures vertical variations due to the above magnetic field 
geometry and b{s) = 4s(rc — s)/r^ takes into account additional magnetic field variations 
with s. The profiles f{x) and b{s) are represented in Figure]^ and the assumed toroidal field 
Figure]^ (left). 


m 


The Elsasser number 




controls the strength of the magnetic pumping, where Bq is the maximum of the amplitude 
of the large scale magnetic field, p is the fluid density and the rotation rate of the Earth. 

Note also that we have approximated the local vorticity by neglecting the contri¬ 

bution of the radial derivative of ^(s, (p). This approximation allows us to conveniently satisfy 
the mass-conservation by adding a contribution to the azimuthal flow: 


1 


= AVb ism^rnis) b{s) j^f{z/H{s)) 


( 11 ) 


Our magnetic pumping flow correction defined by equations Q and 0 enhances the 
helicity of the flow while satisfying the mass conservation as well as the impenetrable boundary 
condition at the core-mantle boundary. It is also noteworthy that the magnetic pumping does 
not alter the surface flow. 

Although we could have parameterized the small correction due to buoyancy forces in a 
similar manner, Sreenivasan and Jones 2011 mention that it has a much smaller effect on 


helicity. Therefore, and to keep our study simple, we ignore the small deformation of columns 
induced by buoyancy effects. 


4 Results 

We have computed kinematic dynamos with a combination of the two parameters A and Rm. 
Using only the mean flow u+ averaged in time and a purely columnar velocity field v given by 
equations!^ we found no dynamo for any tested value of Rm up to Rm = 2600. This is not 
unexpected as Schaeffer and Cardin [2006 also did not find dynamo action with stationary 
columnar flows. 

In addition, adding the time dependence prescribed by our principal component analysis 
did not help, and no dynamo was found up to Rm = 2600. Moreover, we checked that the 
average growth rate is exactly the same as for the mean flow only. Willis and Gubbins 2004| 


showed that adding low frequency perturbations to a mean flow did not change the growth 
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Figure 3: Radial (left) and vertical (right) profiles of our parametrized magnetic pumping. 
They were chosen to capture the dependence with s and z corresponding to a simple field of 
dipolar symmetry [see figure lb of 


Sreenivasan and Jones 2011 


rate, while higher frequency often enhanced dynamo action. Since our principal components 
have periods larger than 60 years, it does not exclude that lower period motions can participate 
in the dynamo process. 

Nevertheless, as an auspicious result, adding magnetic pumping (eq. and 11) to our 
core flows allowed them to produce growing magnetic fields: the larger the Elsasser number 
A, the lower the critical magnetic Reynolds number, down to Rnic ~ 250. Below A = 0.5, we 
could not find a dynamo for Rm up to 1700. These computations are summarized in table 
and figure showing a regime diagram for dynamo action with magnetic pumping. We note 
that the time-dependence has still no effect on the average growth rate, even in the presence 
of magnetic pumping. 

Since our flow v is symmetric with respect to the equatorial plane, the magnetic field eigen 
modes of the kinematic dynamo problem also have a definite symmetry (either symmetric or 
anti-symmetric). It turns out that the fastest growing magnetic field is anti-symmetric with 
a large dipolar component (see figure]^, which is consistent with the assumption used to 
model the magnetic pumping (see section [3.3.3 ). 

Because the growing magnetic field is likely different from the simple assumed toroidal field 
responsible for the magnetic pumping in our framework, the resulting kinematic dynamos are 
not a priori dynamically self-consistent. As an a posteriori consistency check, we can compare 
the growing axisymmetric toroidal magnetic field to the one we have used to compute the 
magnetic pumping in equations and 11 As shown in Fig. they are in fair agreement, 
although more complexity can be seen in the growing magnetic field. In particular they are 
both equatorially antisymmetric. 
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Figure 4: Growing radial magnetic field shown at the top of the core (Aitoff projection). It 
has been rescaled to match the assumed Elsasser number A = 0.5 (left), A = 0.9 (middle), 
A = 2 (right). For each A, the top row is the largest Rm whereas the bottom row is the 
smallest Rm leading to a growing field (see tablefor precise values). The units are /rT 





Figure 5: Growth rate as a function of the magnetic Reynolds number Rm and the Elsasser 
number A controlling the magnetic pumping. Black squares are failed dynamos (for which 
the magnetic field decays) while red circles are dynamos (exhibiting growing magnetic field). 
The contour lines are interpolated using the shown data points. The growth rate has been 
normalized using the magnetic diffusion time. The data is given in table 
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A 

Rm 

Rm* 

7 

BrmsI Bdip 

0.15 

1727 

600 

-7.6 


0.25 

1733 

602 

-0.4 


0.3 

1303 

453 

-1.28 


0.3 

1502 

522 

1.15 

18.1 

0.3 

1737 

604 

4.44 

17.2 

0.4 

915 

318 

-1.5 


0.4 

1107 

385 

1.5 

13.6 

0.4 

1308 

455 

5.3 

12.9 

0.4 

1744 

606 

14.9 

14.7 

0.5 

800 

279 

- 0.01 


0.5 

700 

244 

-1.7 


0.5 

1000 

348 

4 

11.2 

0.5 

1400 

488 

14.5 

12.1 

0.5 

1751 

610 

27 

14.0 

0.7 

503 

176 

-2 


0.7 

596 

209 

-0.08 


0.7 

706 

247 

2.55 

8.4 

0.7 

1059 

371 

12.7 

10.1 

0.9 

400 

142 

-1 


0.9 

489 

173 

0.74 

6.5 

0.9 

649 

230 

4.66 

7.4 

0.9 

978 

346 

14 

9.2 

1.2 

252 

91.2 

-1.9 


1.2 

360 

130 

2.3 

5.5 

1.2 

495 

179 

6.4 

6.4 

1.5 

201 

74.7 

-2.9 


1.5 

222 

82.5 

- 1.6 


1.5 

274 

102 

1.4 

5.1 

1.5 

411 

153 

8 

6.2 

2 

222 

80.7 

- 1.1 


2 

252 

92 

0.9 

5.2 

2 

454 

165 

11 

6.9 


Table 1: Kinematic dynamo runs using the mean core flow and our magnetic pumping param¬ 
eterization. The magnetic Reynolds number is computed using the maximum velocity (Rm) 
or the rms velocity {Rm*). 7 is the growth rate in magnetic diffusion time units, computed 
from the temporal evolution of the magnetic energy Esit) as 7 = {2Eb)~^ dtEs- For the 
growing magnetic field cases (7 > 0 ), the amplitude ratio between the root-mean-square field 
in the bulk and the dipole at the surface is given by B^ms/Bdip- 
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Figure 6: Left: the simple magnetic field used in eqs. and [m to model the magnetic 
pumping. Right: axisymmetric part of the growing magnetic field for Rm = 978 and A = 0.9. 
The colormap shows the azimuthal (toroidal) component and the black contours show the 
meridional (poloidal) field lines. The small gray circle marks the size of Earth’s inner-core, 
although we have none in our study. 

5 Discussion 


While Guervilly et ah 2012 used a forcing at the surface to produce a dynamic bulk flow 


compatible with the observed zonal jets of giant planets, here we decided to entirely prescribe 
the flow, as the core flow is more complex and time dependent. 

We have observed growing magnetic fields in columnar flows with magnetic pumping, for 
magnetic Reynolds number Rm = Vorcji] as low as 252 (for Elsasser number A = 2). If we use 
the root-mean-square velocity in the bulk instead of the maximum velocity Vq, it translates 
to Rm* = 92. This rather low value of Rm* indicates that magnetic pumping in columnar 
flows is a rather efficient mechanism for dynamo action. 

Our flow is a model of the large scale flow in the Earth’s core where, because the operating 
dynamo is saturated, the Lorentz force backreacts on the flow to prevent further exponential 
growth of the magnetic field. Flows from such saturated dynamos have been found to still 


be efficient kinematic dynamos Cattaneo and Tobias, 2009, Tilgner and Brandenburg 


If our flow is mainly responsible for the saturated geodynamo, then according to Cattaneo 


2008 


and Tobias [2009 it is also a kinematic dynamo for a growing magnetic field not everywhere 
proportional to the geomagnetic field and not backreacting on the flow (a passive field). Our 
kinematic dynamo calculations thus exhibit the property of velocity fields arising from a 
saturated dynamo to lead to exponential growth of a passive magnetic field. 

Furthermore, in a kinematic dynamo, the obtained held is an eigenmode, growing or 
decaying at the same rate everywhere, effectively enslaving the time variation of the surface 
held to the time variation in the bulk. It may be worth emphasizing that because of vanishing 
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Figure 7: Growing radial magnetic field (colormap) and streamlines of the surface mean flow 
for Rm = 978 and A = 0.9. Top: Aitoff projection at the core surface centered on the pacific; 
bottom; north-pole view projected onto the equatorial plane (the Greenwich meridian is on 
the right, as in Fig. [^. The thickness of the streamlines is proportional to the velocity. The 
dashed circle marks the boundary of the Earth’s inner-core. 


14 













Figure 8: Anti-symmetric part of the radial magnetic field of gufml, averaged from 1840 to 
1990 (Hammer projection). Units are in /rT. 


radial velocities at the core surface, the magnetic field there is connected to the bulk field only 
through magnetic diffusion. As a result, in a regime of exponentially growing field, the dipole 
moment is also expected to grow. In contrast, the present geomagnetic field is expectedly 
in a saturated regime, where Lorentz forces play an important role back-reacting on the 
flow, although diffusion has probably also a non negligible effect. In particular, the observed 
rapid dipole moment decrease has been ascribed to either growth of reverse geomagnetic flux 
patches in the Southern Hemisphere [e.g. Gubbins et al., 2006 , or advection of normal flux to 

2006]. Hence, particular 


the equator and of reversed flux to the poles [e.g. Olson and Amit 


features of the observed SV are much dependent on the particular geodynamo regime and 
cannot be considered here. In particular, we cannot compare the growth rate obtained for 
our passive field with the current decay rate of the geomagnetic (saturated) held strength. 
In fact, the present decay of the zonal dipole magnetic held energy which is being observed 
for the last 170 years, could in principle ht into the time variability observed in saturated 
geodynamos [see e.g. Christensen, 2011 hg. 6]. 


Cattaneo and Tobias "[2009 have also shown that, in spite of their differences, the growing 


passive held and the self-consistent saturated held have similar spatial structures. It may thus 
be interesting to compare the spatial structure of the growing magnetic held that emerges 
from our kinematic dynamos to the Earth’s internal magnetic held observed nowadays. To 
ease this comparison, the mean anti-symmetric magnetic held of gufml is shown in hgure 
[^ The morphology of the growing held does not vary too much within the portion of the 
parameter space we have explored (see Fig. |^. It is always of dipolar symmetry, consistent 
with the simple toroidal held assumed to compute the magnetic pumping (eq[^. As a general 
trend, larger Rm have more pronounced small-scale features, and large A lead to fewer but 
stronger intense hux patches, and also lower amplitude reverse hux patches near the equator. 
A noticeable feature is the strong hux patch located right under the Bering Strait, which is 
associated with a large vortex in the mean how, as can be seen in hgure [^ 

We want to emphasize that a perfect match with the Earth’s magnetic held is not expected. 
Indeed, our mean how computed over 150 years is unlikely to be an accurate estimate of 
the mean how over a magnetic diffusion time (about 10® yrs). Hence, the strong magnetic 
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flux patch under the Bering Strait, clearly associated with an intense vortex in the mean 
flow as seen in figure and having no corresponding feature in the present Earth’s field, 
may indicate that this very vortex is not part of the long-term mean flow. However, the 
presence of symmetric flux patches in the Earth’s field not far from this location (with northern 
hemisphere signatures under North America and Eastern Russia, fig. |^, may also suggest 
that the whirl structure has suffered some displacement around its present position, over very 
long time periods. The fact that this vortex is also present in the first EOE (see figure 
supports our hypothesis. 

It is also satisfying that, when rescaling the magnetic field to match the imposed A, the 
radial field at the top of the core has the same order of magnitude as the Earth’s field. Eur- 
thermore, the amplitude ratio between the root-mean-square field in the bulk and the dipolar 
field at the surface is around ~ 10 (see table [^, in agreement with numerical geodynamo 
models [Christensen 2011 and torsional oscillations of ~ 6 yrs periodicity retrieved from 
geomagnetic field models Gillet et al. 2010| . 


Another striking result is the footprint of the inner-core in the magnetic field generated 
deep inside the sphere, in spite of no inner core being present in our numerical computations. 
This display of the inner core presence is manifest in the change of sign of the magnetic field 
near the (virtual) tangent cylinder (see figure]^. This is related to the visible dichotomy in 
the prescribed flows (see figures and [^ . 

The mechanism of dynamo generation that is considered in this study involves helicity. 
Here, an enhancement of helicity is prescribed in the bulk, which propagates no footprint 
whatsoever up to the core surface. It is, accordingly, contrasting with the scenario consid¬ 
ered by Amit and Olson 2004j where the in-phase occurrence of vorticity and flow shear in 
the volume is assumed to manifest on the core surface in the form of correlated horizontal 
divergence and radial vorticity. 

Eor A = 0.25 and below, no dynamos were found for Rm < 1700, as shown in figure 
It appears that below some magnetic field strength, dynamo action due to the magnetic- 
pumping is lost, meaning that a strong magnetic field is needed for this mechanism to work. 
This is known as a subcritical behaviour Sreenivasan and Jones 2011 Morin and Dormy 


2009] . Our findings thus suggest that today, a subcritical dynamo may actually produce the 
magnetic field of our planet. If this is the case, a sudden change in our planet’s field may 
occur when in the distant future the magnetic field bulk intensity drops below some threshold. 

The magnetic-pumping model used in this study is admittedly crude. A more realis¬ 
tic model, which would dynamically adjust the magnetic-pumping to the actual magnetic 
field, would allow more realistic kinematic dynamos, with dynamic adjustment of the pump¬ 
ing flow to the growing magnetic field. It may also be possible to improve the quasi- 
geostrophic dynamos of Schaeffer and Cardin 2006] and Guervilly [2010 towards non-linear 
quasi-geostrophic geodynamo models. 


6 Conclusion 

We have shown that magnetic pumping in columnar core flows, leads to dynamo action. The 
distinctive feature of our study is to use a prescribed velocity field and its time evolution, 
derived from current geomagnetic field observations. This expectedly gives them an Earth- 
like nature, exhibiting features with specific geographical location and vorticity distribution 
that are most relevant to test their dynamo action. Note also that our kinematic dynamo 
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approach allows us to reach realistic values of the electrical conductivity. 

Furthermore, we show that relaxing the strict axial invariance with a magnetic-pumping is 
enough for large scale core flows to produce a dipolar field that resembles the one of the Earth, 
without the contribution of small scales. As an additional consistency check, the magnetic 
field in the bulk is about 10 times larger than at the surface (see table [^, in agreement with 
the study of Gillet et al. 2010 . Magnetic pumping is arguably an important mechanism for 


the geodynamo. 

In our computations, the magnetic field generation takes place in the bulk of the core 
where the flow is perturbated, and the surface magnetic field is connected to the bulk field 
only by diffusion. In the Earth’s core, in a saturated but fluctuating regime, it is still unclear 
how much magnetic diffusion can contribute to the SV. Studies bridging numerical dynamo 
simulations and geomagnetic field inversions may help to shed more light on the subject, by 
providing estimates of the radial diffusion at the core surface. 

Eor about fifty years, time-dependent geomagnetic field models have been used to infer 
the geometry and intensity of large scale flows in the core responsible for secular variation. 
Are these flow features expected to be present at diffusive timescales, much larger than the 
time window where we can constrain them with observations? The answer depends also on 
their ability for dynamo action. Our present results, together with previous studies proposing 
a possible dynamical mechanism for the main mean flow features such as the large eccentric 
jet [see Aubert et al., 2013 , suggest indeed that they could be present for a very long time. 
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